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ABSTRACT 


A  high-resolution,  multi-level,  primitive  equation  ocean  model  is  used  to  examine  the 
response  of  the  coastal  region  from  22.5°  S  to  35°  S  of  the  Chile  Current  System  to  both 
equatorward  and  climatological  wind  forcing.  The  results  from  both  types  of  forcing 
show  that  an  equatorward  surface  current,  a  poleward  undercurrent,  upwelling,  meanders, 
filaments  and  eddies  develop  in  response  to  the  predominant  equatorward  wind  forcing. 
When  climatological  wind  forcing  is  used,  an  offshore  branch  of  the  equatorward  surface 
current  is  also  generated.  These  features  are  consistent  with  available  observations  of  the 
Chile  Current  System.  The  model  results  support  the  hypothesis  that  wind  forcing  is  an 
important  mechanism  for  generating  currents,  eddies  and  filaments  in  the  Chile  eastern 
boundary  current  system  and  in  other  eastern  boundary  current  regions  which  have 
predominantly  equatorward  wind  forcing. 
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I.  INTRODUCTION 


The  eastern  boundary  current  (EBC)  system  off  the  west  coast  of  Chile  is  part  of 
the  Peru  Humboldt  coastal  upwelling  system,  which  extends  from  5°  S  to  50°  S.  This 
system  is  influenced  predominantly  by  equatorward  winds  throughout  the  year  (Figure  1). 
The  direction,  and  to  a  large  degree  the  magnitude  of  the  winds,  are  the  result  of  a  semi¬ 
permanent  subtropical  high  pressure  system  that  is  similar  in  nature  and  behavior  to  its 
counterpart  in  the  Northern  Hemisphere,  the  North  Pacific  Subtropical  High,  as  described 
in  Nelson  (1977).  The  South  Pacific  Subtropical  High  migrates  meridionally  with  the 
seasons,  reaching  its  northernmost  extent  during  the  austral  winter  (Figures  lb,  lc),  with 
its  center  located  roughly  at  25°  S.  As  the  seasons  progress,  the  high  pressure  system 
moves  southward  until  its  center  reaches  about  32.5°  S  during  the  austral  summer 
(Figures  le,  If). 

Because  of  this  migration,  maximum  values  of  wind  stress  (Figure  1)  vary 
temporally  at  given  locations.  North  of  Chile,  for  example,  the  maximum  wind  stress  is 
seen  during  the  austral  winter  (Figures  lb,  lc),  presumably  because  the  South  Pacific 
Subtropical  High  reaches  its  northernmost  migration  at  that  time,  and  the  resulting  winds 
take  on  more  of  an  alongshore  characteristic.  The  middle  of  the  Chilean  coastline,  in 
contrast,  does  not  experience  maximum  wind  stress  values  until  late  spring  and  summer 
(Figures  Id,  le). 

The  oceanic  regime  adjacent  to  the  coast  of  Chile  is  comprised  of  several  major 
features  (see  Figure  2),  including  an  equatorward  surface  current  called  the  Chile-Peru 
Current  (CPC)  and  a  poleward  undercurrent  known  as  both  the  Gunther  Current  and  the 
Peru-Chile  Undercurrent  (PCU).  The  climatological  CPC  is  similar  to  other  EBCs  in  that 
it  is  a  sluggish  (Order  (10  cm/s)),  wide  and  generally  shallow  flow  moving  equatorward 
along  a  west  coast  (Wooster  and  Reid,  1963).  Mass  transports  have  been  estimated  to  be 
~19  x  106  m3/s  (Wooster  and  Reid,  1963).  The  primary  source  of  the  CPC  is,  as 
expected,  the  West  Wind  Drift,  because  the  South  American  continent  extends  far  enough 
south  to  intercept  at  least  part  of  the  flow,  and  divert  it  towards  the  north.  The  CPC 
consists  of  two  equatorward  currents  (Tchernia,  1980;  Gunther,  1936).  The  outer  current 


1 


is  usually  referred  to  as  the  oceanic  CPC,  which  has  a  maximum  depth  of  ~  700  m 
(Tchernia  1980).  The  inner  surface  current  is  known  as  the  coastal  CPC,  which  extends 
to  ~  200  m  depth  (Tchernia,  1980). 

Until  recently,  the  PCU  had  been  relatively  unexplored.  Although  Gunther  (1936) 
discussed  the  undercurrent  that  would  eventually  carry  his  name,  very  little  additional 
research  was  conducted  on  the  PCU  until  1961,  when  Wooster  and  Gilmartin  (1961) 
measured  it  using  parachute  drogues.  Their  article  sparked  a  new  interest  in  the  region, 
and  soon  several  expeditions  set  out  to  study  the  waters  adjacent  to  the  South  American 
coast.  Between  1960  and  1982  the  area  was  explored  by  twelve  different  expeditions 
(Fonseca.  1989).  Results  of  the  expeditions  revealed  that  the  PCU  is  present  all  year,  has 
its  core  at  -150-200  m  depth,  is  principally  located  between  the  coast  and  200  km 
offshore,  has  average  speeds  on  the  order  of  10  cm/s,  and  has  a  mass  transport  of —10  x 
106  m3/s  (Fonseca,  1989). 

Recent  observations  have  shown  that,  superimposed  on  the  broad,  climatological 
mean  flow  in  this  and  other  EBC  regions,  such  as  the  California  Current  System,  are 
highly  energetic,  mesoscale  features  such  as  meanders,  eddies  and  filaments  (Mooers  and 
Robinson,  1984).  These  features  have  been  observed  during  periods  of  predominantly 
equatorward  winds,  which  are  favorable  for  upwelling.  These  observations  provide 
evidence  for  wind  forcing  as  a  possible  important  mechanism  for  the  formation  of 
currents,  meanders,  eddies,  and  filaments  in  EBC  regions. 

The  Chile  Current  System  offers  a  unique  opportunity  to  isolate  the  effects  of 
predominantly  equatorward  winds  in  an  EBC  system.  Since  the  coastline  is  relatively 
straight  (see  Figure  3),  there  should  be  no  significant  effects  of  wind  forcing  with 
irregular  coastline  geometry  (such  as  capes  and  bays).  Also,  since  the  topography 
abruptly  drops  off  the  coast  to  the  Peru-Chile  Trench,  there  should  be  no  major  effects  of 
wind  forcing  with  topography. 

In  contrast,  the  other  EBC  regions  that  have  predominantly  equatorward  winds, 
such  as  the  central  and  southern  regions  of  the  California  Current  System,  have  irregular 
coastline  geometry  and  topography.  Also,  in  other  EBC  regions  where  the  coastlines  are 
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relatively  straight,  such  as  in  the  northern  regions  of  the  California  Current  System  and  of 
the  Canary  Current  System,  the  wind  stress  is  predominantly  equatorward  only  during  the 
upwelling  season,  from  around  April  to  September  (During  the  rest  of  the  year,  the  winds 
reverse  and  become  predominantly  poleward.).  In  the  EBC  region  off  Western  Australia, 
the  winds  are  predominantly  equatorward  and  the  coastline  is  relatively  straight;  however, 
the  thermal  forcing  is  so  strong  in  this  region  that  the  predominantly  observed  surface 
current  (the  Leeuwin  Current)  opposes  the  winds  (e.g..  Batteen  et  al.  1992b). 

There  have  been  numerous  wind  forcing  models  used  over  the  past  few  decades  to 
model  EBCs,  particularly  the  California  Current  System.  Early  work  included  that  of 
Pedlosky  (1974)  and  Philander  and  Yoon  (1982),  who  used  steady  wind  stress  and 
transient  wind  forcing,  respectively.  Carton  (1984)  and  Carton  and  Philander  (1984) 
investigated  the  response  of  reduced  gravity  models  to  realistic  coastal  winds.  Another 
series  of  experiments  were  done  by  McCreary  et  al.  (1987)  utilizing  a  linear  model  with 
both  transient  and  steady  wind  forcing  in  the  California  Current  System.  In  all  of  these 
models,  relatively  weak  currents  (5-10  cm/s)  were  generated  but  no  eddies  or  filaments 
developed. 

Batteen  et  al.  (1989)  used  a  primitive  equation,  multi-level  model  with 
biharmonic  (rather  than  Laplacian,  as  used  by  the  previous  authors)  heat  and  momentum 
diffusion  which  had  both  steady  and  meridionally  varying  wind  forcing.  In  both  cases, 
when  the  baroclinic/barotropic  shear  between  the  surface  current  and  undercurrent 
became  strong  enough,  meanders  and  filaments  developed.  Although  this  was  the  first 
EBC  model  to  simulate  meanders  and  filaments  in  the  California  Current  System,  no 
eddies  were  generated  by  the  end  of  the  90-day  experiments,  presumably  due  to  the 
relatively  short  simulation  time. 

The  objective  of  this  study  is  to  extend  the  work  of  Batteen  et  al.  (1989)  by 
investigating  the  effects  that  both  equatorward  and  climatological  wind  forcing  have  on 
the  generation  of  not  only  currents,  meanders,  and  filaments,  but  also  eddies,  in  the  Chile 
Current  System.  A  slightly  modified  version  of  the  high-resolution,  multi-level,  primitive 
equation  model  of  Batteen  et  al.  (1989)  will  be  used  to  simulate  the  effects  of  wind  stress 
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in  this  region.  Longer  simulation  times  than  90  days  will  be  used  to  allow  the  eddies  to 
be  generated  and  subsequently  analyzed. 

The  organization  of  this  study  is  as  follows:  Chapter  II  describes  the  numerical 
model,  the  forcing,  the  specific  initial  and  experimental  conditions,  and  the  energy 
analysis  technique.  This  technique  is  used  to  investigate  the  dynamical  reasons  for  the 
generation  of  eddies  in  the  Chile  Current  System.  Chapter  III  includes  an  analysis  of  the 
results  of  the  modeling  experiments,  along  with  a  comparison  of  model  results  with 
available  observations,  while  Chapter  IV  contains  a  summary. 
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II.  THE  MODEL 


A.  MODEL  DESCRIPTION  AND  INITIAL  CONDITIONS 

To  investigate  the  role  of  wind  forcing  on  the  generation  of  currents,  eddies,  and 
filaments  in  the  Chile  EBC  region,  the  wind  stress  fields,  discussed  below,  were  used  to 
specify  the  wind  forcing  for  a  high-resolution,  multi-level,  primitive  equation  (PE)  model 
of  a  baroclinic  ocean  on  a  (3-plane.  The  model  is  based  on  the  hydrostatic,  Boussinesq, 
and  rigid  lid  approximations.  The  governing  equations  are  defined  in  Batteen  et  al. 
(1989).  For  the  finite  differencing,  a  space-staggered  B-scheme  (Arakawa  and  Lamb, 
1977)  is  used  in  the  horizontal.  Batteen  and  Han  (1981)  have  shown  that  this  scheme  is 
appropriate  when  the  grid  spacing  is  approximately  on  the  same  order  as,  or  less  than,  the 
Rossby  radius  of  deformation,  which  meets  the  criteria  of  this  study.  The  horizontal  grid 
spacing  is  10  km  alongshore  and  9  km  cross-shore,  while  the  internal  Rossby  radius  of 
deformation  is  ~30  km.  In  the  vertical,  the  10  layers  are  separated  by  constant  z  levels  at 
depths  of  13,  46,  98,  182,  316,  529,  870,  1416,  2283,  and  3656  m.  Consistent  with  Haney 
(1974),  this  spacing  scheme  concentrates  more  on  the  upper,  dynamically  active  part  of 
the  ocean,  above  the  thermocline. 

The  model  domain  (Figure  3)  extends  off  the  west  coast  of  Chile,  from  -22.5°  S 
to  35°  S  (1280  km  alongshore),  and  from  -72.5°  W  to  82.5°  W  (1152  km  cross-shore). 
The  eastern  boundary  of  the  model  domain  is  closed,  and  has  both  the  tangential  and 
normal  components  of  velocity  set  to  zero.  To  isolate  the  role  of  wind  forcing  in  the 
generation  of  eddies,  the  eastern  boundary  is  modeled  as  a  straight,  vertical  wall 
extending  from  the  surface  to  4500  m  depth.  These  assumptions  are  justified  since  the 
Chilean  coastline  is  relatively  straight  and  smooth  (see  Figures  2  and  3),  the  continental 
shelf  is  relatively  narrow,  and  the  continental  slope  is  steep,  dropping  from  sea  level  to 
over  2000  m  depth  in  a  distance  of  less  than  25  km  from  the  coast  (Ziegler  et  al.  1957). 
The  northern,  southern  and  western  borders  are  open  boundaries  which  use  a  modified 
version  of  the  radiation  boundary  conditions  of  Camerlengo  and  O’Brien  (1980).  Some 
spatial  smoothing  is  applied  in  the  vicinity  of  the  open  boundaries. 
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The  model  uses  biharmonic  lateral  heat  and  momentum  diffusion  with  the  same 
choice  of  coefficients  (i.e.,  2.0  x  1017cmV)  as  in  Batteen  et  al.  (1989).  Holland  (1978) 
showed  that  the  highly  scale-selective  biharmonic  diffusion  acts  predominantly  on 
submesoscales,  while  Holland  and  Batteen  (1986)  found  that  baroclinic  mesoscale 
processes  can  be  damped  by  Laplacian  lateral  heat  diffusion.  As  a  result,  the  use  of 
biharmonic  lateral  diffusion  should  allow  mesoscale  eddy  generation  via  barotropic 
(horizontal  shear)  and/or  baroclinic  (vertical  shear)  instability  mechanisms.  As  in  Batteen 
et  al.  (1989),  weak  (0.5  cmV1)  vertical  eddy  viscosities  and  conductivities  are  used. 
Bottom  stress  is  parameterized  by  a  simplified  quadratic  drag  law  (Weatherly,  1972),  as 
in  Batteen  et  al.  (1989). 

The  method  of  solution  is  straightforward  with  the  rigid  lid  and  flat  bottom 
assumptions  because  the  vertically  integrated  horizontal  velocity  is  subsequently 
nondivergent.  The  vertical  mean  flow  can  be  described  by  a  stream  function  which  can 
be  predicted  from  the  vorticity  equation,  while  the  vertical  shear  currents  can  be  predicted 
after  the  vertical  mean  flow  is  subtracted  from  the  original  equations.  The  other 
variables,  i.e.,  temperature,  density,  vertical  velocity  and  pressure,  can  be  explicitly 
obtained  from  the  thermodynamic  energy  equation,  equation  of  state,  continuity  equation, 
and  hydrostatic  relation,  respectively.  Since  density  is  primarily  a  function  of 
temperature  (e.g.,  see  Wooster  and  Reid,  1963)  and  there  are  no  major  salinity  sources  or 
sinks  (such  as  major  rivers)  in  the  region  being  modeled,  effects  of  slinity  changes  on 
density  are  neglected. 

An  annual  climatological  temperature  profile  (Figure  4),  centered  at  -28°  S 
(corresponding  to  the  middle  of  the  model  domain),  based  on  Levitus  (1982),  was  used  to 
represent  the  initial  temperature  conditions  for  the  model  domain.  An  exponential  curve 
was  then  used  to  approximate  this.  The  form  of  the  equation  used  was: 

T(z)  =  T,i  +  A  Te:/h 

The  approximation  uses  a  length  scale  of  450  m,  and  assumes  TB=2°  C  to  be  the 
temperature  at  great  depth.  AT=15°  C  is  the  change  of  temperature  between  the  bottom 
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and  top  of  the  profile.  Figure  4  shows  that  the  resulting  temperature  profile  fits  the 
annual  climatological  temperature  profile.  The  exponential  profile  was  used  to  provide 
initial  values  at  each  of  the  ten  levels  in  the  model. 

B.  SPECIFIC  EXPERIMENTAL  DESIGN 

1.  Experiment  1:  Constant  Equatorward  Winds 

In  the  first  experiment,  a  constant  equatorward  wind  stress  of  1  dyne/cm~  (~8 jO 
cm/s)  is  applied  to  all  grid  points  within  the  domain,  except  for  a  small  band  at  the 
northern  and  southern  ends,  where  no  wind  forcing  is  applied.  McCreary  (1981)  and 
Batteen  et  al.  (1989)  demonstrated  that  this  technique  allowed  for  the  propagation  of 
coastal  Kelvin  waves,  which  in  turn  allowed  for  the  creation  of  an  undercurrent.  Because 
the  winds  are  predominantly  equatorward  all  year  off  Chile,  this  experiment  should 
isolate  the  effects  of  upwelling  favorable  winds  on  the  Chile  Current  System. 

2.  Experiment  2:  Variable  Climatological  Winds 

In  Experiment  2,  climatological  winds  are  used  to  simulate  the  temporal  and 
spatial  wind  variations  that  have  been  observed  off  the  coast  of  Chile.  The  winds  used  to 
force  the  model  are  from  the  European  Center  for  Medium-Range  Weather  Forecasts 
(ECMWF)  (Trenberth  et  al.  1990)  and  have  been  provided  by  the  National  Center  for 
Atmospheric  Research.  Winds  were  extracted  from  world  wide  mean  monthly  wind 
velocities  that  were  provided  at  2.5  degree  spacing  for  the  years  1980  to  1989  along  with 
120  month  climatology  for  the  area  of  study.  This  data  was  linearly  interpolated  in  space 
and  then  in  time  to  provide  daily  forcing  at  all  of  the  grid  points. 

An  example  of  the  wind  fields  used  can  be  seen  in  Figure  5,  which  shows  the 
annual  migration  of  the  South  Pacific  Subtropical  High  from  the  south  in  the  austral 
summer  (Figure  5d),  to  the  north  in  the  austral  winter  (Figure  5b),  the  effects  of  which  are 
discussed  in  the  next  chapter.  These  wind  field  patterns  are  consistent  with  those  of 
Bakun  and  Nelson  (1991). 

To  isolate  the  effects  of  wind  forcing  in  this  process-oriented  study,  the  net  heat 
flux  at  the  sea  surface  should  initially  be  zero.  In  this  way,  any  heat  flux  that  is 
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discernible  will  have  been  generated  by  an  increase  or  decrease  of  sea  surface  temperature 
resulting  from  wind  forcing  effects  (Batteen  et  al.  1989).  To  accomplish  this,  an  initial 
air  temperature  is  chosen  that  forces  the  net  flux  of  longwave  radiation,  sensible  heat,  and 
latent  heat  to  zero.  This  air  temperature  is  then  used  in  the  model  for  both  experiments. 
Any  subsequent  surface  heat  flux  forcing  is  therefore  a  secondary  effect  of  the  changes  to 
sea  surface  temperature  due  to  the  wind  forcing. 

C.  ENERGY  ANALYSIS  TECHNIQUE 

The  energy  technique  used  is  the  same  as  that  used  and  described  in  Batteen  el  al. 
(1992b),  and  is  based  on  that  of  Han  (1975)  and  Semtner  and  Mintz  (1977).  This  analysis 
is  done  to  gain  a  better  understanding  of  the  types  of  energy  transfer  during  unstable  flow 
in  the  Chile  Current  System.  A  brief  summary  of  the  method  follows. 

Kinetic  energy  is  calculated  for  the  horizontal  components.  After  quasi-steady 
state  is  reached  where  the  total  kinetic  energy  is  nearly  constant,  mean  and  eddy  kinetic 
energies  are  calculated  using  the  averaged  sum  of  squared  mean  and  eddy  horizontal 
fields,  respectively.  Next,  the  available  potential  energy  is  calculated  and  used  to 
determine  when  a  quasi-steady  state  is  reached  and  when  statistics  should  be  collected. 
Then,  both  mean  and  eddy  available  potential  energies  are  computed.  The  barotropic  and 
baroclinic  energy  transfers,  defined  in  Batteen  et  al.  (1992b),  are  used  to  argue  for  the 
type  of  instability  mechanism  (e.g.,  barotropic,  baroclinic,  or  mixed)  which  leads  to  the 
initial  eddy  generation  in  each  experiment. 
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III.  RESULTS  OF  THE  EXPERIMENTS 


Experiments  1  and  2  study  the  role  of  wind  forcing  in  the  formation  of  currents 
and  eddies  in  the  Chile  Current  System.  Experiment  1  uses  steady  equatorward  winds, 
while  Experiment  2  uses  climatological  average  (1980-1989)  winds. 

A.  EXPERIMENT  1  -  CONSTANT  EQUATORWARD  WINDS 

Experiment  1  was  forced  with  a  constant  equatorward  wind  stress  (1  dyne/cm2), 
which  was  uniform  in  both  the  alongshore  and  cross-shore  directions.  As  expected,  the 
equatorward  wind  forcing  resulted  in  an  equatorward  surface  current  along  most  of  the 
coast  (Figure  6a).  The  current’s  velocities  ranged  from  ~10  to  30  cm/s  throughout  the 
duration  of  the  experiment.  The  coastal  jet  axis  was  within  ~50  km  of  the  coast  and 
extended  to  -150  m  depth  (Figure  7). 

A  poleward  undercurrent  with  a  core  velocity  of -8  cm/s  at  ~  100  to  150  m  depth 
is  also  seen  in  Figure  7  below  the  surface  current.  The  undercurrent  extends  along  the 
entire  coast  and  is  confined  to  ~50  km  of  the  coast.  Typical  velocities  for  the 
undercurrent  ranged  from  ~  5  to  15  cm/s  throughout  the  duration  of  the  experiment. 

Along  with  the  equatorward  surface  current,  upwelling  also  occurred,  as  seen  in 
Figure  8a.  Meanders  in  the  equatorward  surface  current  formed  by  model  day  45  and  are 
evident  in  both  the  velocity  and  temperature  plots  (Figures  6b  and  8b,  respectively).  As 
the  meanders  intensify,  cold  upwelling  filaments  develop  along  the  coast  and 
subsequently  extend  farther  offshore  (Compare  day  45  in  Figure  8b  with  day  75  in  Figure 
8c  to  see  the  time  evolution  ). 

More  detailed  analysis  was  performed  to  determine  the  type  of  instability 
mechanism  that  could  generate  the  meander  features.  Barotropic  instability  can  result 
from  horizontal  shear,  while  baroclinic  instability  can  result  from  vertical  shear  in  the 
current.  As  Figure  7  shows,  there  is  considerable  horizontal  and  vertical  shear  between 
the  equatorward  surface  current  and  the  poleward  undercurrent.  As  a  result,  both  types  of 
instabilities  (mixed),  can  be  present  simultaneously.  Energy  transfer  calculations  which 
consist  of  barotropic  (mean  kinetic  energy  to  eddy  kinetic  energy)  and  baroclinic  (eddy 
potential  energy  to  eddy  kinetic  energy)  components  were  performed  for  the  time  period 
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(i.e.,  days  39-60)  during  which  the  meanders  developed.  The  results  for  the  instability 
analyses  (Figure  9)  show  that  both  barotropic  (Figure  9a)  and  baroclinic  (Figure  9b) 
instabilities  were  present  in  the  coastal  region  of  the  domain. 

By  model  day  105,  the  meanders  began  to  form  cold  core,  cyclonic  eddies  and 
warm  core,  anticyclonic  eddies,  as  seen  in  Figures  6c  and  8d,  due  to  vertical  and 
horizontal  shear  instabilities  between  the  equatorward  jet  and  the  poleward  undercurrent. 
The  eddies  initially  appeared  to  be  about  100  km  in  diameter  and  extended  about  100-200 
km  off  the  coast.  Subsequently,  the  eddies  propagated  farther  offshore  with  speeds  of  ~5- 
10  km/day,  consistent  with  Rossby  wave  propagation  speeds. 

At  the  end  of  the  model  year,  the  coastal  jet,  undercurrent  (not  shown),  meanders 
and  eddies  (which  have  propagated  further  offshore)  are  still  present  (Figure  6d). 
Upwelling  and  filaments  are  also  still  discernible  (Figure  8e).  These  features  are  due  to 
the  constant  equatorward  wind  forcing,  which  provides  a  continual  source  for  the 
generation  of  the  currents  and  the  subsequent  development  (via  barotropic/baroclinic 
instability  mechanisms)  of  meanders,  filaments,  and  eddies.  Longer  experimental  runs  of 
three  more  years  (not  shown)  show  that  these  features  continue  to  be  generated  and 
maintained  in  the  Chile  Current  System. 

B.  EXPERIMENT  2  -  VARIABLE  CLIMATOLOGICAL  WINDS 

Experiment  2  was  forced  with  the  climatological  average  (1980-1989)  winds 
(Figure  5).  An  equatorward  surface  current  was  observed  all  along  the  coast  by  model 
day  30  (Figure  10a)  due  to  the  prevailing  equatorward  winds.  The  current’s  velocities 
ranged  from  -10  to  30  cin/s  throughout  the  experiment.  The  coastal  jet  axis  was  within 
~50  km  of  the  coast  and  extended  to  -250  m  depth  (Figure  1 1).  These  results  are  similar 
to  Experiment  1,  as  expected,  since  the  winds  are  still  predominantly  equatorward. 

Unlike  Experiment  1,  an  equatorward  surface  current  also  develops  -200  km 
offshore  between  days  30  and  45  in  the  poleward  end  of  the  model  domain  (Figure  10b). 
The  core  of  the  current  extends  to  a  depth  of  -250  m  depth  (Figure  11).  From  days  90  to 
150,  it  extends  along  the  entire  coast  (e.g.,  see  Figure  10c).  From  days  165  to  210,  it  is 
present  only  in  the  equatorward  end  of  the  domain  (e.g.,  see  Figure  lOd).  It  subsequently 
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redevelops  in  the  poleward  end  of  the  domain  by  day  225  (Figuie  lOe).  The  location  ot 
this  feature  appears  to  migrate  with  the  South  Pacific  Subtropical  High,  which  migrates 
from  the  south  in  the  austral  summer  (Figure  5d)  to  the  north  in  the  austral  winter  (Figure 

5b). 

In  addition  to  the  coastal  and  offshore  branches  of  the  surface  equatorward 
currents,  a  poleward  undercurrent  with  a  core  velocity  of -3  cm/s  at  ~1 50  m  depth  is  also 
seen  in  Figure  11  below  the  coastal  equatorward  current.  As  in  Experiment  1,  the 
undercurrent  extends  along  the  entire  coast  and  is  confined  to  ~50-75  km  of  the  coast. 
Typical  velocities  for  the  undercurrent  range  from  3  to  15  cm/s  throughout  the  duration  of 
the  experiment. 

Along  with  the  equatorward  surface  currents,  upwelling  also  occurs,  as  seen  in  the 
surface  temperature  field  (Figure  12a).  The  coastal  location  of  colder  temperatures  due  to 
upwelling  appears  to  migrate  from  poleward  (Figure  12a)  to  equatorward  (Figure  12d), 
possibly  following  the  migration  of  the  South  Pacific  Subtropical  High. 

Meanders  in  the  equatorward  surface  current  formed  by  model  day  135  in  the 
equatorward  end  of  the  model  domain,  as  seen  in  the  velocity  and  temperature  plots 
(Figures  10c  and  12b,  respectively).  The  results  for  the  instability  analysis  (Figure  13) 
for  the  time  period  (i.e.,  days  120-135)  during  which  the  meanders  developed,  show  that, 
although  both  types  of  instabilities  are  present  in  the  coastal  region  of  the  model  domain, 
barotropic  instability  (Figure  13a)  is  dominant  over  baroclinic  (Figure  13b)  instability  in 

the  equatorward  end  of  the  model  domain. 

By  model  day  1 80,  the  meanders  begin  to  form  cold  core,  cyclonic  eddies  as  seen 
in  the  temperature  and  velocity  plots  (Figures  lOd  and  12c,  respectively),  due  to  vertical 
and  horizontal  shear  instabilities  between  the  equatorward  jet  and  the  poleward 
undercurrent.  By  model  day  225,  the  meanders  also  begin  to  form  warm,  anticyclonic 
eddies,  as  seen  in  the  temperature  and  velocity  plots  (Figure  12d  and  lOe,  respectively). 
Both  the  cyclonic  and  anticyclonic  eddies  initially  appear  to  be  about  100  km  in  diameter 
and  extend  about  200  km  off  the  coast.  As  in  Experiment  1,  the  eddies  propagate  farther 
offshore  with  speeds  of  -5-10  km/day.  Unlike  Experiment  1,  as  the  eddies  move 
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offshore,  they  become  embedded  in  the  oceanic  branch  of  the  surface  equatorward 
current,  which  veers  further  offshore  with  time  (see  Figures  lOe  and  lOf,  for  example). 

At  the  end  of  the  model  year  these  features  are  still  discernible.  Upwelling  and 
filaments  are  still  present  (Figure  12e).  Longer  experimental  runs  of  three  more  years 
(not  shown)  show  that  these  features  continue  to  be  generated  and  maintained  in  the  Chile 
Current  System. 

The  coastal  branch  of  the  equatorward  surface  current,  the  undercurrent, 
upwelling,  meanders,  filaments,  and  eddies  are  due  to  the  predominantly  equatorward 
wind  forcing.  These  features  are  similar  to  Experiment  1.  Unlike  Experiment  1,  an 
offshore  branch  of  the  equatorward  surface  current  is  also  generated.  This  feature 
originates  in  the  vicinity  of  the  South  Pacific  Subtropical  High  and  its  location  migrates 
with  the  location  of  the  South  Pacific  Subtropical  High. 

C.  COMPARISON  OF  MODEL  RESULTS  WITH  OBSERVATIONS 

A  comparison  of  model  results  with  available  observations  was  carried  out  to  see 
if  model  simulations  of  the  currents  and  other  features  were  consistent  with  the  observed 
data  for  the  Chile  Current  System.  Table  1  highlights  the  primary  features  in  the  Chile 
Current  System.  It  also  lists  the  specific  characteristics  (if  available  observations  exist  for 
them)  for  each  feature  along  with  the  corresponding  model  result  for  each  experiment. 

Table  1  shows  that  both  experiments  successfully  simulate  most  of  the  observed 
features  of  the  Chile  Current  System,  i.e.,  the  coastal  equatorward  surface  current  (the 
coastal  CPC),  the  coastal  poleward  undercurrent  (the  PCU),  upwelling,  filaments, 
meanders,  and  eddies.  Consistent  with  observations,  the  experiments  also  successfully 
simulate  the  magnitudes  of  the  average  velocity,  core  depth,  and  mass  transport  for  the 
coastal  CPC,  and  the  average  velocity,  core  depth,  offshore  distance,  and  mass  transport 
for  the  PCU.  Only  Experiment  2  successfully  simulates  the  offshore  equatorward  surface 
current  (the  offshore  CPC).  In  addition.  Experiment  2  also  successfully  simulates  the 
magnitudes  of  the  average  velocity,  core  depth,  and  offshore  distance  for  the  offshore 
CPC. 
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IV.  SUMMARY 


This  study  used  a  high  resolution,  multi-level,  primitive  equation  model  to 
investigate  the  role  of  wind  forcing  in  the  formation  of  currents  and  eddies  in  the  Chile 
Current  System.  Experiment  1  used  steady  equatorward  winds,  while  Experiment  2  used 
climatological  average  (1980-1989)  winds  obtained  from  the  European  Center  for 
Medium-Range  Weather  Forecasts  (ECMWF). 

The  model  results  of  both  experiments  showed  the  development  of  a  coastal 
equatorward  surface  current  and  coastal  poleward  undercurrent.  Due  to 
barotropic/baroclinic  instability,  the  coastal  currents  became  unstable,  resulting  in  the 
generation  of  meanders  near  the  coast.  Upwelling  also  occurred  and,  as  the  meanders 
intensified,  cold  upwelling  filaments  developed  along  the  coast  and  subsequently 
extended  farther  offshore.  In  time,  the  meanders  formed  both  cold  core,  cyclonic  eddies 
and  warm  core,  anticyclonic  eddies,  which  subsequently  propagated  farther  offshore  with 
speeds  consistent  with  Rossby  wave  propagation  speeds.  By  the  end  of  the  model  year, 
all  of  these  features  were  still  present,  as  expected,  since  the  predominantly  equatorward 
winds  used  in  both  experiments  provide  a  continual  source  for  the  generation  of  the 
currents  and  the  subsequent  development  of  upwelling,  meanders,  filaments,  and  eddies. 

Only  the  model  results  of  Experiment  2  showed  the  development  of  an  offshore 
equatorward  surface  current.  The  location  of  this  feature  appeared  to  migrate  with  the 
South  Pacific  Subtropical  High,  which  migrates  from  the  south  in  the  austral  summer  to 
the  north  in  the  austral  winter. 

A  comparison  of  model  results  with  available  observations  showed  that  both 
experiments  successfully  simulated  most  of  the  observed  features  of  the  Chile  Current 
System,  i.e.,  the  coastal  CPC,  the  PCU,  upwelling,  filaments,  meanders,  and  eddies.  Only 
Experiment  2  successfully  simulated  the  offshore  CPC,  which  is  another  observed  feature 
of  the  Chile  Current  System. 

The  results  from  these  experiments  strongly  support  the  hypothesis  that  wind 
forcing  can  be  a  significant  generation  mechanism  for  currents,  upwelling,  filaments, 
meanders,  and  eddies  in  the  Chile  Current  System.  The  results  of  this  study  may  also  be 
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applicable  to  other  EBC  regions  which  have  predominantly  equatorward  wind  forcing, 
such  as  the  central  and  southern  regions  of  the  California  Current  System. 


14 


Figure  1 .  Wind  stress  fields  for  the  time  period  (a)  April-May,  (b)  June- 
July,  (c)  August-September,  (d)  October-November,  (e)  December- 
January,  and  (f)  February-March.  Vectors  are  scaled  in  units  of  dyne/cm  . 
Ellipses  show  areas  of  large  standard  errors  (after  Bakun  and  Nelson, 
1991). 
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Figure  2.  Generalized  circulation  schematic  for  the  Chile  Current  System 
(after  Codispoti  et  al.  1989):  The  broad  equatorward  Chile-Peru  Current 
(CPC)  (solid  lines)  separates  into  two  branches:  the  coastal  CPC  and  the 
oceanic  CPC.  The  coastal  CPC  overlies  the  poleward  undercurrent,  known 
as  the  Peru-Chile  Undercurrent  (dashed  line). 


Latitude  (degrees) 

i  i 


Figure  3.  Model  Domain:  Area  of  study  is  a  1152  km  (cross-shore)  by 
1280  km  (alongshore)  box  off  the  coast  of  Chile. 
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TEMPERATURE  PROFILE 

ANNUAL 


Temperature  (degrees  Celsius  ) 

Figure  4.  Temperature  profile  (i.e.,  temperature  in  degrees  Celsius  versus 
depth  in  meters)  used  in  the  model:  A  smooth  exponential  curve  (solid 
line)  was  used  to  approximate  the  climatological  temperature  profile  (solid 
line  with  small  circles)  computed  from  Levitus  (1982). 
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Figure  5.  Climatological  (1980  -  1989)  ECMWF  winds  in  m/s:  (a)  day  90, 
(b)  day  180,  (c)  day  270,  (d)  day  360.  The  latitudinal  grid  point  1  (129) 
corresponds  to  35°  S  (22.5°  S).  The  longitudinal  grid  point  1  (129) 
corresponds  to  82.5°  W  (72.5°  W).  Maximum  wind  vector  is  20  m/s. 


19 


Distance  off  shore  (km) 


Distance  off  shore  (km) 


Velocity  at  model  day  30.0 


Velocity  at  model  day 


g.egge^2 
MAXIMUM  ‘Vrtros 


Distance  off  shore  (km) 


Distance  off  shore  (km) 


Js'-Y, 


r_  ».»- 
A  =* 

"  V  o 

V  o 

\  Vf  .c 
/  i  '  o> 

, .  Yi  *  j  6W  °- 

_v . < 


/  c 


/  c 


r~  v  W 


£r££T  Y  t 


✓  /  »  / 


Ya' 


Velocity  at  model  day  105.0 


Velocity  at  model  day  360.0 


(c)  (d) 

Figure  6.  Experiment  1 :  Surface  velocity  vectors  in  the  coastal  region  at 
days  (a)  30,  (b)  45,  (c)  105,  and  (d)  360.  To  avoid  clutter,  velocity  vectors 
are  plotted  at  every  third  grid  point  in  both  the  cross-shore  and  alongshore 
directions,  and  velocities  less  than  5  cm/s  are  not  plotted.  The  numbers 
associated  with  the  highs  in  Figures  6c  and  6d  correspond  to  temperature, 
not  velocity.  Maximum  current  vector  is  80  cm/s. 
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Figure  7.  Experiment  1:  Cross-shore  section  at  y~640  km  (~29°S)  of  the 
meridional  component  of  velocity  at  day  30.  Solid  lines  indicate 
equatorward  flow,  while  dashed  lines  indicate  poleward  flow.  The  contour 
interval  is  2  cm/s.  The  maximum  equatorward  velocity  is  ~20  cm/s  at  the 
surface,  ~50  km  off  the  coast. 
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Figure  8.  Experiment  1:  Surface  temperature  in  the  coastal  region  at  days 
(a)  30,  (b)  45,  (c)  75,  (d)  105,  and  (e)  360.  The  contour  interval  is  0.5°  C. 
The  temperature  decreases  toward  the  coast.  Note  the  existence  of 
upwelling  and  meanders. 
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Figure  9.  Experiment  1 :  Energy  transfers  in  the  coastal  region  of  (a)  mean 
to  eddy  kinetic  energy  (i.e.,  barotropic  energy  transfer)  and  (b)  eddy 
potential  to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer)  for  model 
days  39  to  60.  The  contour  interval  is  10  ergs  cm"3 s’1. 
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Figure  10.  Experiment  2:  Surface  velocity  vectors  in  the  coastal  region  at 
days  (a)  30,  (b)  45,  (c)  135,  (d)  180,  (e)  225,  and  (f)  360.  To  avoid  clutter, 
velocity  vectors  are  plotted  at  every  third  grid  point  in  both  the  cross-shore 
and  alongshore  directions,  and  velocities  less  than  5  cm/s  are  not  plotted. 
The  numbers  associated  with  the  highs  in  Figures  lOd  and  lOf  correspond 
to  temperature,  not  velocity.  Maximum  current  vector  is  80  cm/s. 
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Figure  11.  Experiment  2:  Cross-shore  section  at  y~  640  km  (-29°  S)  of 
the  meridional  component  of  velocity  at  day  30.  Solid  lines  indicate 
equatorward  flow,  while  dashed  lines  indicate  poleward  flow.  The  contour 
interval  is  2  cm/s.  The  maximum  equatorward  velocity  is  ~20  cm/s  at  the 
surface.  ~50  km  off  the  coast. 
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Figure  12.  Experiment  2:  Surface  temperature  at  days  (a)  45,  (b)  135,  (c) 
180,  (d)  225,  and  (e)  360.  The  contour  interval  is  0.5°  C.  The  temperature 
decreases  toward  the  coast.  Note  the  existence  of  upwelling  and 
meanders. 
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Figure  13.  Experiment  2:  Energy  transfers  in  the  coastal  region  of  (a) 
mean  to  eddy  kinetic  energy  (i.e.,  barotropic  energy  transfer)  and  (b)  eddy 
potential  to  eddy  kinetic  energy  (i.e.,  baroclinic  energy  transfer)  for  model 
days  120  to  135.  The  contour  interval  is  10  ergs  cm'3 s'1. 
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j  Features: 

Observation  j  Experiment  1 

Experiment  2 

A  Coastal  Equatorward  Surface  Current  (coastal  Peru-Chile  Current  (CPC)): 

1  Average  Velocity  (cm/s) 

>10(1) 

10-30 

10-30 

2  Core  Depth  (m) 

200  (2) 

150 

250 

3  Mass  Transport  (10fa  m'Vs) 

19(1) 

20 

20 

J  B  Offshore  Equatorward  Surface  Current  (offshore  CPC):  [ 

1  Average  Velocity  (cm/s) 

10-15  (2) 

not  simulated 

10 

2  Core  Depth  (m) 

<700  (2) 

not  simulated 

250 

|  3  Offshore  Distance  (km) 

200  (2) 

not  simulated 

200 

C  Coastal  Poleward  Undercurrent  (Peru-Chile  Undercurrent  (PCU)):  j 

1  Average  Velocity  (cm/s) 

10(3) 

5-15 

3-15 

2  Core  Depth  (m) 

150-200  (3) 

100-150 

150 

3  Offshore  Distance  (km) 

<200  (3) 

50 

75-100 

4  Mass  Transport  (10bmd/s) 

10(3) 

10 

10 

D  Upwelling 

observed 

simulated 

simulated 

E  Filaments 

observed 

simulated 

simulated 

F  Meanders 

observed 

simulated 

simulated 

G  Eddies 

observed 

simulated 

simulated 

Table  1.  Comparison  Of  Model  Results  With  Available  Observations 
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